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We define a large class of continuous time multifractal random measures and processes with 
arbitrary log-infinitely divisible exact or asymptotic scaling law. These processes generalize within a 
unified framework both the recently defined log-normal Multifractal Random Walk (MRW) [jj], [| and 
the log-Poisson "product of cynlindrical pulses" J3 . Our construction is based on some "continuous 
stochastic multiplication" (as introduced in [Q) from coarse to fine scales that can be seen as a 
continuous interpolation of discrete multiplicative cascades. We prove the stochastic convergence 
of the defined processes and study their main statistical properties. The question of genericity 
(universality) of limit multifractal processes is addressed within this new framework. We finally 
provide a method for numerical simulations and discuss some specific examples. 
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I. INTRODUCTION 

Multifractal processes are now widely used models 
in many areas including nonlinear physics, geophysics 
or econophysics. They are used to account for scale- 
invariance properties of some observed data. Our pur- 
pose in this article is to introduce a wide class of random 
measures and stochastic processes with stationary incre- 
ments that possess exact multifractal scaling without any 
preferred scale ratio. This construction generalizes (and 
unifies) the recently introduced log-normal "Multifractal 
Random Walk" model jl], ||] and the log-Poisson com- 
pounded "Multifractal Product of Cynlindrical Pulses" 
]3|. Technical mathematical proofs are reported in the 
companion paper || . 

In the multifractal framework, scale-invariance prop- 
erties of a Id stochastic process X(t) are generally char- 
acterized by the exponents £ g which govern the power 
law scaling of the absolute moments of the "fluctuation" 
SiX(t) of X(t) at any scale I up to a scale T, i.e., 

m(q,l) = K q l c «,Vl<T, (1) 

where m(q, I) is defined as the expectation (from now 
on, the symbol E{.) will always refer to the mathematical 
expectation) : 

m(q,l) = E[\6 l X(t)\«]. (2) 

The so-called "fluctuation" SiX(t) can be defined in 
various ways, the most commonly used definition being 

5iX(t)=X(t + l)-X(t). (3) 



When the exponent is linear in q, i.e., 3£T, £ 9 = qH, 
the process is referred to as a monojractal process. Let us 
note that the so-called self-similar processes (e.g., frac- 
tionnal Brownian motions, Levy walks) are a particular 
case of monofractal processes. On the contrary, if C, q is a 
non-linear function of q it is referred to as a multifractal 
process (or as a process displaying multiscaling or inter- 
mittency). 

The scale-invariance property (]]]) can be qualified as a 
continuous scale-invariance in the sense that the relation 
is a strict equality for the continuum of scales < I < T. 
Alternative "weaker" forms of scale-invariance have been 
widely used. All of them assume at least the asymptotic 
scale-invariance relation 

m(q,l) ~ K q f<>, l^0 + . (4) 

The discrete scale-invariance property adds strict equal- 
ity for discrete scale values l n (with l n — > + , when 
n — * +oo) 

m(q, l n ) = K q l%>, n — > +oo. (5) 

The image implicitely associated with a multifractal 
process is a random multiplicative cascade from coarse to 
fine scales. Such cascading processes can be explicitely 
defined in very different ways ranging from the original 
construction proposed by Mandelbrot in its early work 
H to recent wavelet-based variants ||. Though these 
different constructions do not lead to the same objects, 
they all aim at building a stochastic process X(t) which 
fluctuation 6\i X (t) at a scale XI (where I is an arbitrary 
scale smaller than the large scale T and A < 1) is obtained 
from its fluctuation SiX(t) at the larger scale I through 
the simple "cascading" rule 
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where W\ is a positive random variable independent of 
X and which law depends only on A. 

In the multiplicative construction by Mandelbrot || 
(which log-normal variant was originally introduced by 
Yaglom [0), a positive multifractal measure M(dt) is 
built, i.e., X(t) = M([0,t]) is an increasing process and 
the fluctuation 6iX(t) at scale / is defined by Eq. (||). 
In the more recent wavelet-based constructions X(t) is 
a process, not necessarily increasing, and the fluctua- 
tion SiX(t) is defined as the wavelet coefficient at scale I 
and time t ||. Such cascade processes have been exten- 
sively used for modeling scale-invariance properties as ex- 
pressed in Eq. (jl]). However, apart from self-similar pro- 
cesses which are monofractal processes with continuous 
dilation- invariance properties, all cascade processes in- 
volve some arbitrary discrete scale ratio and consequently 
only have discrete scale-invariance properties (Eq. (||)). 
Indeed, they rely on a "coarse to fine" approach consist- 
ing first in fixing the fluctuations at the large scale T and 
then, using Eq. (^) iteratively, deriving fluctuations at 
smaller and smaller scales. However, the (t, I) time-scale 
half-plane is very constrained : one cannot choose freely 
SiX(t) at all scales and times. Thus, the cascade is gen- 
erally not built using the whole time-scale half-plane but 
only using a sparse (e.g., dyadic) grid in this half-plane. 
Consequently an arbitrary scale ratio (e.g., A = 1/2) is 
introduced in the construction. The continuous scale- 
invariance property (|IJ) is thus broken and replaced by 
the weaker discrete scale-invariance property (^j) with 
l n = \ n T. Moreover, as a consequence, the fluctuations 
6iX(t) in these approaches are no longer stationary. 

Let us point out that despite the potential interest 
of stochastic processes with stationary fluctuations and 
continuous invariance-scaling properties, until recently 
jl], ||, [|, explicit constructions possessing such proper- 
ties were lacking. Let us assume that one can build a 
continuous cascade process X(t) satisfying Eq. (^) with 
a continuous dilation parameter A. As first pointed 
out by Novikov if such a construction is possible, 
a simple transitivity argument shows that ln W\ must 
have an infinitely divisible law: its characteristic func- 
tion is of the form G\(q) — \ F ( q \ One can then easily 
prove that m(q, XI) = m{q,l)\ Fty ~ tq \ and consequently 
C, q = F{—iq), Continuous cascade statistics and log- 
infinite divisibility have been the subject of many works 
with applications in various domains ranging from turbu- 
lence to geophysics §, [jl], HI |J, |J, [[! H, || . In 
the case W\ = \ H is deterministic, one gets £ 9 = qH and 
therefore X(t) is a monofractal process. This is the case 
of the so-called self-similar processes. The simplest non- 
linear (i.e., multifractal) case is the so-called log-normal 
cascade that corresponds to a log-normal law for W\ and 
thus to a parabolic £ q spectrum. Other well known log- 
infinitely divisible models often used in the context of 
fully developed turbulence are log-Poisson |18 and log- 
Levy Q models. 

From our knowledge, among all the attempts to build 
multifractal processes with continuous scale invariance 



properties and stationary fluctuations, only the recent 
works by Bacry et al. jl], |j and Barral and Mandelbrot 
H refer to a precise mathematical construction. Bacry 
et al. have built the so-called Multifractal Random Walk 
(MRW) processes as limit processes based on discrete 
time random walks with stochastic log-normal variance. 
Independently, Barral and Mandelbrot || have proposed 
a new class of stationary multifractal measures. Their 
construction is based on distributing, in a half-plane, 
Poisson points associated with independent identically 
distributed (i.i.d.) random weights and then taking a 
product of these weights over conical domains. In the 
same time, Schmitt and Marsan considered an extension 
of discrete cascades to a continuous scale framework Q] 
and introduced infinitely divisible stochastic intregrals 
over cone-like structures similar to those of Barral and 
Mandelbrot. However they did not considered any con- 
tinuous time limit of their construction nor studied its 
scaling properties as defined by Eq. (]!]). 

In this paper, we propose a model that is based on the 
stochastic approach developed in [Q and that unifies all 
previous constructions within a single framework. Start- 
ing from original discrete multiplicative cascades, we will 
use a cone-like construction as in Q and || in order to 
get rid of discrete scale ratios and to consider any log- 
infinitely divisible multifractal statistics. We will show 
that this allows us to build a very large class of multi- 
fractal measures and processes (including original MRW 
El H and Barral-Mandelbrot multifractal measures jj|) 
for which both long range correlations and multiscaling 
properties can be controlled very easily. 

The paper is organized as follows. In section II, we re- 
view the discrete multiplicative cascades in order to nat- 
urally introduce the notion of stochastic integral over a 
cone-like structure in some "time-scale" half-plane. We 
then define a class of log-infinitely divisible stationary 
Multifractal Random Measures (MRM) which statisti- 
cal properties are studied in section III. In section IV, 
we define the log-infinitely divisible Multifractal Random 
Walks. We show how their scaling properties can be in- 
ferred from the associated MRM. In section V we address 
some questions related to numerical simulations and pro- 
vide explicit examples. In section VI, we discuss some 
links of the present work with previous connected ap- 
proaches. Conclusions and prospects for future research 
are reported in section VII. 



II. FROM DISCRETE MULTIPLICATIVE 
MEASURES TO MULTIFRACTAL RANDOM 
MEASURES 

In this section, we provide some heuristics about 
how one can build a positive stationary stochastic mea- 
sure M(dt) with continuous scale-invariance properties, 
i.e., such that the associated increasing process X(t) = 
M([0,i]) satisfies Eq. (0) with a continuous dilation pa- 
rameter A. 



A. Discrete multiplicative cascades 

For the sake of illustration, let us start with simple 
discrete multiplicative cascades. In the original construc- 
tion, as proposed by Mandelbrot, one builds the measure 
M(dt) as the limit of a sequence of stochastic measures 
M in , indexed by a discrete scale parameter l n = TX n 
(we choose A = 1/2). The measure Mi n , at the step n 
of the construction, has a uniform density on successive 
dyadic intervals of size /„. The idea is to build the se- 
quence Mi n so that it statisfies (when n is varying) a 
scale-invariance property. The n-th measure is obtained 
from the (n— l)-th measure by multiplication with a pos- 
itive random process W which law does not depend on 
n. One can naturally index the dyadic intervals along the 
dyadic tree using a kneading sequence {si . . . s n } where 
Si = (resp. Si — 1) if, at "depth" i, the interval is on 
the left (resp. on the right) boundary of its parent inter- 
val. Thus, for instance, one gets the following dyadic 
intervals: Jo = [0,2- x ], h = [2" 1 , 1], 7 00 = [0,2" 2 ], 
Zoi = [2~ 2 ,2 -1 ] and so on... With these notations, the 
multiplicative rule reads (see Fig. 1(a)): 

M ; „(7 S0 ... S J = W^.^M^I.,...^) (7) 

where I So ... Srl is of size T2~ n and is one of the two "sons" 
of the interval J So ... Sn _ 1 and W ao ,., Sn are i.i.d. Since the 
construction is invariant with respect to a rescaling by 
a factor 2, the limit measure M(dt) will then satisfy the 
same scale-invariance property and will be multifractal 
in the discrete sense ofEq. (§). There is a huge math- 
ematical literature devoted to the study of such a con- 
struction and we refer the reader to Refs. ^2| |2^] 
for rigourous results about the existence, regularity and 
statistical properties of Mandelbrot cascades. In physics 
or other applied sciences, as recalled in the introduction, 
the previous construction (and many of its variants) is 
considered as the paradigm for multifractal objects and 
has been often used as a reference model in order to re- 
produce observed multiscaling. But because of its lack of 
continuous scale invariance and translation invariance, 
such models cannot be fully statisfactory in many con- 
texts were the considered phenomena possess some de- 
gree of stationarity and do not display any preferred scale 
ratio. 
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FIG. 1: Defining a continuous cascade interpolating Mandel- 
brot cascades, (a) Mandelbrot cascade: One starts from the 
coarse scale T and constructs the sequence Mi n that is re- 
solved at scale T2~" by iterating a multiplicative rule along 
the dyadic tree. At each construction step, the random multi- 
plicative factor W ao ...a n can be represented as the exponential 
of the integral of an infinitely divisible noise over a square 
domain indicated by hatched sets, (b) Continuous cascade 
conical domains obtained as linear smoothing of the hatched 
domains in (a). Such cones are involved in the definition of 
the process ui 2 -3 (t) 



B. Revisiting discrete cascades 

In order to generalize Eq. (^) to a continuous frame- 
work one can try to perform the limit A — > 1 in the dis- 
crete construction as in Ref. Q| . Another way to proceed 
is to represent a Mandelbrot cascade as a discretization 
of an underlying continuous construction. 

For that purpose, we suppose that the random weights 
W So ,.. Sn in the cascade are log-infinitely divisible. For a 
continuous cascade, this choice can be motivated as fol- 
lows: Let us suppose that the large scale density, Mt (dt) 



is equal (or proportional) to the Lebesgue measure dt. 
We would like to define iteratively the densities Mi(t), 
for all I <T. To go from resolution V to resolution I < I' , 
we thus write 

Mi(t) l =Wi /v (t)M v {t), VI < I', (8) 

where Wi/i'(t) is a positive stationary discrete random 
process independent of M[{t). Let us define 

w,(t)=lnW I/T (i). (9) 
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Thus, one gets, 

Mi(t)dt = e" l(t) M T (t)dt = e^dt. (10) 

By iterating Eq. (^J), we see that Wj/x(i) can be obtained 
as a "continuous product" of positive i.i.d. random vari- 
ables. Consequently u>i(t) can be written as the sum of 
an arbitrary number of i.i.d. random variables. This is 
precisely the definition of an infinitely divisible random 
variable j24|. In order to go from discrete to continu- 
ous cascades, it is therefore natural to assume that, in 
the discrete situation, w;„(so ■ ■ ■ s n ) — In (]X_o W So ... Si ) 
is infinitely divisible. A "simple" way of building such 
an infinitely divisible process 0Ji n is to represent it by a 
stochastic integral of an infinitely divisible stochastic 2D 
measure P(dt, dl) over a domain Ai n (so ■ ■ • 



Lj ln (so...s n )=P(At n {sa...s n ))= / P(dt,dl) 

J Ai n (s ...s„) 

(11) 

The stochastic measure P(dt, dl) is uniformly (with re- 
spect to a measure fx(dt, dl)) distributed on the time-scale 
half-plane S + = {(t, I), i e R, le M+}. 

Let us recall that, by definition 1 , if P(dt,dl) is a 
stochastic infinitely divisible measure uniformly (with re- 
spect to a measure fi(dt, dl)) distributed on S + , then, for 
any two ^i-mesurable sets A and A' such that (J,(A) = 
fJ,(A'), then P(A) and P(A') are identically distributed 
random variables which characteristic function is nothing 
but 



(12) 



where f(p) depends only on a centering parameter m 
and the so-called canonical Levy measure v(dx) which is 
associated with P. The general shape of ip is described 
by the celebrated Levy-Khintchine formula pj, [2B] : 



ip(p) = imp 



e tpx — 1 — ip sin x 



v(dx). (13) 



with J_ v(dx)/x 2 < 00 and J v(dx)/x 2 < 00 for all 

y>o. 

The sets Ai n (so . . . s n ) associated with each of the 
2 n values Wj„=t2-« ( s o ■ ■ ■ s n ) in Mandelbrot construction 
can be chosen naturally as the union of all similar squares 
I S0 ... Sk x [T2- fc ,T2-( fe - 1 )], n < k < 1: 

Ai n (s ... sn) = ^ =0 Is ...s k x [T2- fe , T2-< fe - 1 )] . (14) 

The sets „4 2 - 3 (000) and ^2-3 (010) are indicated as 
hatched domains in Fig. 1(a). Since the W / s ... s „'s are 
i.i.d., we want to choose the measure fi(dt,dl) such that 



the measure of each square I So ...s„ x [T2~ n , T2~ (n ~ 1 )] is 
a constant. The natural measure to choose is fi(dt,dl) = 
dtdl/l 2 . It is the natural measure associated with the 
time-scale plane S + in the sense that it is (left-) invari- 
ant by the translation-dilation group. 

Fixing fi{dt,dl) — dtdl/l 2 one then gets Vn, 



d\x = I r' 2 dl I dt = 1/2 

, ,.„„ x[T2-™,T2-("- 1 >] Jl- 

and thus, from Eq. ([l2|), W SL 
exactly the Mandelbrot construction. 



are i.i.d.: we recover 



Towards multiplicative cascades with 
continuous scale-invariance 



We would like to apply the previous scheme in the case 
the construction is no longer indexed by a discrete scale 
parameter /„ but by a continuous scale parameter I. 

In order to "interpolate" smoothly this construction 
both in time and scale, on can interpolate the previous 
union of similar squares Ai n (sq . . . s n ) (Eq. (jl4)) using 
domains Ai(t) where (t, I) can take any value in S + . In 
order the limit measure M(dt) = lim;^o Mi[dt) to be 
stationary, it is clear that one has to choose the set Ai (t) 
to be "translation-invariant" in the sense that 

(t + T,l')eAi{t + T)^{t,l')eAi(t), Vr. (15) 

One thus just needs to specify the set „ 4;(0). A "natural" 
choice (though, as we will see in section III F , it only leads 
to asymptotic, and not exact, scaling properties) seems 
to be the conical set (see Fig. 1(b)) 

MO) = {(i, l'), l' > l -/(0/2 < * < /GO/2}, (16) 
where f(l) is the function 

f(l) = 



l for I < T 
for / > T, 



(17) 



In case case, we get exactly the construction originally 
proposed, using a densification argument, by Schmitt and 
Marsan in Ref. J5J. Other choic es for the function f(l) 
will be discussed in section [II F . 

Let us remark that the choice of the linear conical 
shape f(l) = I is consistent with the interpretation of 
the parameter / as a scale parameter. Indeed, the value 
of the measure P around some position (to, I) in the half- 
plane S + influences the values of the measure over the 
time interval ]to — 1/2, to + 1/2], i.e., exactly over a time 
scale I. 



III. STATIONARY MULTIFRACTAL RANDOM 
MEASURES 



1 we refer the reader to ref. |26| for a rigourous definition of the so- 
called "infinitely divisible independently scattered random mea- 



A. Denning MRM 

According to the arguments of the previous section, 
we thus propose the following definition for the class of 
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log-infinitely divisible multifractal random measures. Let 
us introduce an infinitely divisible stochastic 2D mea- 
sure P uniformly distributed on the half-plane S + = 
{(t, I), t E R, ! e R + } with respect to the measure 
/i(dt,dl) = dtdl/l 2 and associated with the Levy mea- 
sure v{dx). Let us recall that for any set A C S + , P(A) 
has an infinitely divisible law whose moment generating 
function is 



E 



ip(-ip)fj,(A) 



(18) 



where f{p) is defined by Eq. (|13[). Henceforth, we define 
the real convex cumulant generating function, as, when 
it exists, 

1>(p) = tp(-ip) . (19) 
Let u>i(t) the stationary stochastic process defined by 

u l (t)=P{A l {t)), (20) 
where Ai (t) is the 2D subset of S + defined by 

Mt) = {(t',i'), i'>i, 

where f(l) satisfies: 



/(0/2<f-f </(0/2}, 
(21) 



f(l) 



I for I < T 
T for / > T, 



(22) 



As it will be shown in the following sections, the choice of 
the large scale behavior f(l) = T for I > T is the unique 
one that ensures the convergence of the construction and 
the exact scaling of the limit measure. The cone-like 
domain Ai(t) is indicated as an hatched domain in Fig. 
2. 

We finally define the stochastic positive measure 
Mi(dt) as 



M t {dt) = e Wlit) dt, 



(23) 



meaning that for any Lebesgue measurable set /, one 
has 2 



M,(J) 



(24) 



The MRM M is then obtained as the limit measure 
(the meaning and the existence of this limit will be ad- 
dressed in the next section) 



M(dt) = lim Mi(dt). 



(25) 



Since a simple change in the mean of the stochastic mea- 
sure P would lead to the same measure up to a determin- 
istic multiplicative factor, we will assume, without loss of 
generality, that ip satisfies 



V(i) = o. 



(26) 



2 Since the paths of u)i (t) are continuous to the rieht and limited 
to the left, the integral (|24h is well defined (see |5|) 



T 



u 




-A f (t) 



Timet 



FIG. 2: Conical domain (Eq. (fca)) in the (t, I) half-plane 
involved in the definition of uii (Eq! (boj) ) . 



We can consider some generalizations of the previous 
construction. The first one consists in changing the func- 
tion f(l), i.e., to change the shape Ai the measure P is 
integrated on. The second one consists in changing the 
measure fx(dt, dl), i.e., to change the way the measure P is 
distributed in the half-plane S + . Actually, from Eq. ( p^ ) 
one can easily show that the construction only depends 
on the function n(Ai), consequently, changing the shape 
of /(/) basically amounts changing fi(dt,dl). A simple 
example that illustrates such a freedom is the choice 
n(dt,dl) = dtdl, i.e. fi is nothing but the 2D Lebesgue 
measure. In that case, n{A{) remains unchanged if one 
chooses /(/) = 1/2 in the definition (p2]). The param- 
eter I is no longer a scale but can be interpreted as a 
frequency and thus S + is the time-frequency half-plane. 
Therefore, in the following sections, without loss of gen- 
erality, we choose to fix fi(dt,dl) = dtdl/l 2 (i.e. to work 
within the time -scale half-plane S + ) and we will discuss, 
in section 1IIF , the consequences which rise from other 
choices than (|22[) for the function f{l). 



B. Existence of the limit MRM M(dt) 

In S, we prove, within the framework of positive con- 
tinuous martingales that, almost surely, Mi converges 
to a well defined limit measure when I — > + . More- 
over, we prove that if ^'(1) < 1 then there exists e > 
such that ^(1 + e) < 1 and the moment of order 1 + e, 
E (M ; ([0,£]) 1+£ ), is finite. Then, using the fact that 
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ip(l) = 0, it is straightforward to prove that 

E(M([0,i\) = Urn E(Mi([0,t]) = 1. (27) 

Consequently the limit measure M(dt) is non degener- 
ated (i.e., different from zero). The overall proof is very 
technical and, for this reason, has not been reproduced 
in this paper. 

However, if ip(2) < 1, one can prove (see Appendix A) 
that sup; E(Mi[0, t] 2 ) is bounded and that the sequence 
Mi(dt) converges in the mean square sense. Again, us- 
ing the fact that ^(1) = 0, it follows that M(dt) is non 
degenerated. (Let us note that, as explained in section 
III D, one can prove |5j that assuming tp(2) < 1 basically 
amounts assuming that E (M([0, t]) 2 ) < +oo). 



C. Exact multifractal scaling of M(dt) 



In Ref. ||, Qz(tq,p q ) is computed using a recurrence 
on q. We obtain the following result: 



QKtq.Pq) = e£J= 1 El= 1 a<3.fc)«(t»-* J ) i (33) 



where 



a(j,k) = f(r k ,j) + <p(r k +i,j-i) - ip{r k ,j-x) - <p(r k+hj ), 

(34) 

and 



ELfcP™. fOT k <3 



> = 1 jf"" ' (35) 



for k > j 



Moreover, let us remark that 



j=l k=l \k=l ) 



In order to study the scaling properties of the limit 
measure M(dt), let us establish the scale invariance prop- 
erties of the process u>i(t). 

1. Characteristic function of u)i(t) 

Let q 6 N*, t q = ti,t<2, . ■ ■ ,t q with t\ < t% < . . . < t n 
and p q = p±,p2, ■ . ■ ,p q . The characteristic function of 
the vector {oJi(t m )}i< m < q is defined by 

Q,(tq,p q ) = £(e^^ p ^»). (28) 

Relation (|l8| ) allows us to get an expression for quantities 
like E ^e^™=i a m P(B m )^j wnere {B m } m would be disjoint 

subsets of S + and a m arbitrary numbers. However the 
{Ai(t m )} m in Eq. ( |28| ) have no reason to be disjoint sub- 
sets. We need to find a decomposition of {Ai(t m )} m onto 
disjoint domains. This is naturally done by considering 
the different intersections between these domains. 
Let us define the cone intersection domains as: 

Ai(t,t')=Ai(t)nAi(t'). (29) 

and 

Pl (t)=(i(Ai(0,t)) . (30) 

Using the definition of Ai(t) with the shape of /(/) as 
given by Eq. (p2|), the expression for pi(t) reads: 

r in m + 1 - f \£t<i 

pi(t) = \ ln(f) iiT>t>l , (31) 

[ if t > T 

Notice that pi(t) satisfies the remarkable property (for 
t<TandA<l): 



2. Multifractal properties of M(dt) 

The multifractal scaling properties of the limit MRM 
M (dt) result from the scale invariance property of the 
process u>i{t) that is itself a direct consequence of pre- 
vious exponential expression (|3^) for the characteristic 
function of the process uji and the particular shape of the 
conical domains leading to expression ( |3l| ) . Indeed, using 
these equations together with Eq. ([36]), it can be proven 
that Vn, Vii, ...t n S [0,T] n , Vpi . ..p n G K", one has 

QA/(Atq, Pq ) = A-^^)Q,(tq,p q ) . 

It follows that, for A < 1, the process u>i{t) satisfies, for 
I £ [0, T], the following invariance property: 

{^(tA)}*'^^ + «,(*)}* (37) 

where £l\ is an infinitely divisible random variable (ie., it 
does not depend on t) which is independent of uii (t) and 
which infinitely divisible law is defined by: 

B( e *pn*) = \-v(p) (38) 

We deduce the following scale invariance relationship 
for the sequence of measures Mj([0,t]): 

{M xl ([0,Xt])} t = { f Xt e^' Mdu} t 
Jo 

X{ e^ (An) du} t '= Ae n H / e^ (u) du}t 
Jo Jo 

= Ae n HM ; ([Q,t])} t = Wx{Mi([0,t])} t , (39) 

where W\ = \e n> - is independent of {M\i([0, Xt])} t . 

By taking the limit I — > + , one gets the continuous 
cascade equation for MRM as defined in Eq. (^): 



Pxi(Xt) = Pl (t)- ln(X). (32) 



{M([0,Xt})} t la ^W x {M([0,t])} t , VAe]0,l], (40) 
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where ln(W\) is an infinitely divisible random variable 
independent of {M([0, t])} t . 

The exact multifractal scaling follows immediately: 
Vq e E, Wt < T, , we get 

E[M{%t]f]=K^% (41) 

were the multifractal spectrum and the prefactor read: 



C q = q-4>(<i) 

K q = T- c "E[M([0,T} 



(42) 
(43) 



Let us notice that the moment E (M ([0, f]) 9 ) in Eq. @ 
can be infinite. Conditions for E (M([Q,t]) q ) to be finite 
will be discussed in section III D. 



3. Multi-scaling of correlation functions 



The exact (multi-)scaling law (41) for the absolute mo- 
ments can be easily extended, along the same line, to 
n-points correlation functions. Indeed, let us define the 
n-point correlation function, when it exists, as 

C n {ti, ■ ■ ■ , t n ; ti, . . . , r n ;pi, . . . ,p n ) — 

E [M{[t!,ti + n]) Pl . . . M([t n , tn + r„]) p "] . (44) 

where r k > 0, and t k + r k < tk+i- It is easy to show that 
if t n + t„ — ti < T, then C n is an homogeneous function 
of degree (^ Pk ■ 

C n ({Xt k }; {Xr k }; {p k }) = \^C n ({t k }; {r k }; {p k }) (45) 
with 



P = 

fe=i 



(46) 



This equation extends the scaling law (|41j) for the mo- 
ments to multi-points correlation functions. 



D. Algebraic tails of probability density functions 

Let us note that, if q > 1, 

E[M([0,m = E[(M([0,t/2]) + M([t/2,t})) q ] 
> 2E[M([0,t/2]) q ]. 

Using the multifractal scaling (p]), one gets > 
2(1/2)^ leading to 1 > 2 1- ^. It follows that 



K„ < oo Co > !• 



(47) 



Thus, if q > 1, the q order moment is infinite if £ q < 1. In 
||, we show that the "reverse" implication is also true, 
i.e., 



Let us notice that this infinite moment condition £ g < 1, 
q > 1, is exactly the same as for discrete multiplicative 
cascades established in Refs. [^0[ |2l], Divergence 
of moments for multifractals have also been discussed in 
e.g. § H- 

Therefore, if there exists some value 1 < < oo such 
that 



C 9 , = 1 (49) 



then, 



Prob{M([0,t]) > x} 



when x — > +oo (50) 



The pdf of the measure M([0, t]) is thus heavy tailed with 
a tail exponent that can be, unlike classical a-stable laws, 
arbitrary large. 



E. Examples 

In order to illustrate previous considerations let us con- 
sider some specific examples. 

• Deterministic case: 

The simplest situation is when the Levy v(dx) 
measure is identically zero: This case corresponds 
to the self-similar, monofractal situation where 
''Pip) — P m - The constraint ^(1) = implies 
m = 0, and we thus get the Lebesgue measure. 

• Log-normal MRM: 

The log-normal MRM is obtained when the canon- 
ical measures attributes a finite mass at the ori- 
gin: v(dx) — \ 2 5{x)dx and A 2 > 0. From Eqs. 
( p"3| ) and (|l^), the cumulant generating function is 
that of a normal distribution: ip(p) = pm + X 2 p 2 /2. 
The condition "0(1) = implies the relationship 
m = — A /2. The log-normal spectrum is a 
parabola: 



dn 



(51) 



Let us note that the so-obtained increasing process 
M([0, t]) is the same as the increasing MRW process 
mentioned in the conclusion of [Q. This similarity 
will be further discussed in section VI.A.l. 

• Log-Poisson MRM: 
When there is a finite mass at some finite value 
xo = hi(<5), of intensity A 2 = 7(ln<5) 2 : u{dx) — 
\ 2 8(x — x$). The corresponding distribution is 
Poisson of scale parameter 7 and intensity ln(<5): 
ip(p) = p(m- sin(ln(5))) - 7(1 - S p ). The log- 
Poisson Cq spectrum is therefore exactly the same 
as the one proposed by She and Leveque in their 
cascade model for turbulence fllSll : 



Cq > 1 => K q < OO. 



(48) 



C; p = gm / + 7(1-^) 



(52) 
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where ml is such that £1 = 1. Notice that in the 
limit 5 — > 1~, 7 (In 5) 2 — > A 2 , one recovers the 
log-normal situation. In the original She-Leveque 
model, 7 = 2 and <5 = 2/3 ||. 

• Log-Poisson compound MRM: 
When the canonical measure v(dx) satisfies 
/ v{dx)x~ 2 = C < oo (e.g. u{dx) is concen- 
trated away from the origin), one can see that 
F{dx) = v(dx)x~ 2 jC is a probability measure. In 
that case, 



<p(p) = im'p + C (e 



l)F(dx) 



is exactly the cumulant generating function associ- 
ated with a poisson process with scale C and com- 
pound with the distribution F pifi . Let us now 
consider a random variable W such that InW is 
distributed according to F(dx). It is easy to see 
that J e px F{dx) = E{W P ). It turns out that the 
log-Poisson compound MRM has the following mul- 
tifractal spectrum: 



d pc = qm-C (E(W q ) - 1) 



(53) 



This is exactly the spectrum obtained by Barral 
and Mandelbrot in their construction of "product of 
cylindrical pulses" ||. The similarity between our 
construction and Barral and Mandelbrot construc- 
tion will be further discussed in section VI. A. 2. 

Log-a stable MRM: 

When v(dx) ~ x 1 ~ a dx for < a < 2, one has an 
log a-stable MRM: 

4s 



C = qm-a a \q\ 



(54) 



Such laws have been used in the context of turbu- 
lence and geophysics |l9| . They have been often 
referred to as "universal multifractals" because in- 
stable laws are fixed points of infinitely divisible 
laws under a suitable renormalization procedure. 

Many other families of ( q spectra can be obtained (e.g. 
log-Gamma, log-Hyperbolic,...) for other choices of the 
Levy measure. Let us remark that in the case of a nor- 
mal random variable loi, the function pi{r) is nothing 
but the covariance of the process as introduced in Refs. 
jl], D . This function, that measures the areas of domains 
Ai intersections, is therefore the analog of the covariance 
for general infinitely divisible distributions. The equa- 
tion (|33] ) shows that our construction can be seen as a 
natural extension of Gaussian processes (or multivariate 
Gaussian laws) within the class of infinitely divisible pro- 
cesses (multivariate infinitely divisible laws) in the sense 
that it is completely characterized by a cumulant generat- 
ing functions (p(p) (specifying the mean and the variance 
parameters in the Gaussian case) and a 2-points "covari- 
ance" function p{t) (or covariance matrix). 



F. Asymptotic scaling and universality 

In the last section, we have seen that the choice ( p2[ ) for 
f(l) leads to exact scaling of the moments of the associ- 
ated MRM. In this section we study the scaling behavior 
of the moments for other choices of /(/) 3 . Let us remark 
that /(/) is defined up to a multiplicative constant. This 
just amounts to a choice in the scale of ip(p). 
In the following, f^ e \l) will refer to the "exact scaling" 
choice (^2|) that was made in the previous sections. The 
so-obtained sets in the S + half-plane will be referred to 
as A[ e \t), the u) process as u>\ (t) and the associated 
MRM as M {e \dt). 



1. Large scale perturbation of f(l) 

Rigorous mathematical proofs can be found in 0]. 
Let us first study the case when one builds an MRM using 
a function f(l) which differs from f^(l) only for scales 
larger than a large (fixed) scale L, i.e., 



Id) 



/W(Z) for I < L 
g(l) for I > L, 



(55) 



For the measure Mi(dt) to remain finite, pi(t) must 
be finite and thus we must have f{u)u~ 2 du < 00. 
Therefore, the large scale behavior g(l) must be such 
that, for some e > 0, g(l) — 0(l 1 ~ e ) as I — > +00. An 
example of such large scale modification is the function 
defined by (|l7|), where g — g^ where g^(l) — and 
L = T. This function is the one which was used by 
Barral and Mandelbrot in Ref. ||. Let us first choose 
the particular case g = g^ (L being any strictly posi- 
tive number). The so-obtained sets in the S + half-plane 
will be referred to as Ai (t), the lu process as oj\ ; (i) 
and the associated MRM as M^ s \dt). Since, V t,t', 
(^ (e) (i)\A (s) ( < )) n A[ s) (t') = 0, one has 

u>j e \t) = 4 s \t) + 8 L (t) 

where S L (t) = P(^ e) (t)\^ w (t)) is a process which is 
independent of the process a>, (t) and which does not 
depend on the value of I (as long as I < L). It follows 
that 



£j(M^([0,t]) 9 ) < E fsupe^Wj £ (mW([0,*])«) 
E (m( £ > ([(),£])«) > E Mnf e^Wj E (a/ w ([0, t]) q 



3 We recall that we do not consider other choices for the uniform 
measure fj,(dt, dl) in the half space S + because /(/) and fi(dt, dl) 
are involved in the properties of the limit measure M(dt) only 
through the function pi(t) = fj,(Ai(0,t)). Hence, up to a change 
of variable I' = h(l), one can always set fj,(dt,dl) = dtdl/l~ 2 . 
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Because the process <5l(£) is (right) continuous, 
lim t ^ 0+ sup [0)t] e« a *W = lim^ 0+ mf [0 , 4] = e^°\ 

we get 

E (M«([0,t])«) ^(7,23 (M«([0,i])«) . (56) 

Thanks to the exact scaling of Af( e )([0, i]) (Eq. @),we 
see that, 

E(M^([0,t}y) tZo C q t^ (57) 

where Cq is defined in Eq. (ff2|).Thus, we find that 
M^ s \dt), corresponding to the specific choice g(l) = 
g( a ) (I) = 0, satisfies the asymptotic scale invariance prop- 
erty (Q). If one chooses a different function g(l), using 
exactly the same arguments as above (in which (dt) 
is replaced by M(di), i.e., the measure obtained when 
using the new function g(l)), one can prove B that 

£(M([0 ) t])«)~ o fl l E(jfW([0 1 (])«) , (58) 

and consequently 

E(M([0,t])i) tZo D q C q t^. (59) 

Thus, using any function g(l) in ( ^5| ) (satisfying g(l) = 
0(Z 1 ~ e ) (2 — > +oo)) leads to an MRM measure which 
satisfies the asymptotic scale invariance property (||). 

2. Small scale perturbation of f(l) 

Let us now study the consequences of a small scale 
perturbation of f^- e \l)- Let us suppose that /(/) ~ l a for 
1^0. In that case, 

( l a ~ l if a < 1 
#(0) - < Cst if a > 1 , (60) 
l ~~ >0 [ -ln(Z) if a = 1 

In the first case, a < 1, one can show any moment 
of order 1 + e (e > 0) cannot be bounded. Since 
E(Mi([0,t])) = 1, general martingale arguments can be 
used to prove that M/([0, t]) converges towards the trivial 
zero measure. 

If a > 1, the limit measure is proportional to the 
Lebesgue measure and thus one obtains the trivial 
asymptotic scaling: 

E(M([0,t]Y) -C^ 

Let us now consider the marginal case, f(l) = I + o{l) 
at small scales. From the results of previous section, with 
no loss of generality, as far as the asymptotic scaling is 
concerned, we can suppose that the large scale compo- 
nent is identical to the exact scaling situation: f(l) = T 
for I > T. We can then show that, once again, the MRM 



satisfies a multifractal scaling in the asymptotic sense 
(Eq.©: 

B(M([0,t])«) t ~ o K'/« (61) 

Indeed, let us index all the quantities by the cut-off scale 
I and the integral scale T, i.e., we add an explicit refer- 
ence to the integral scale T: pi^if) is the area of domain 
intersections and Mi t T([0,i\) is the associated MRM. Af- 
ter some little algebra, using the definition of pi^r{t) and 
the fact that /(/) = I + o(Z), one can show that 

PxiMW -^P^kt) . (62) 
Thus, thanks to Eq. ([33]), we have 

Hence, because M l>T ([0, t]) = f* e Ul - T ^du, by taking the 
limit I — > 0, 

A- 1 M O) Ar([0,At])^Mg([0,t]) 
and therefore, from Eqs. ©, © and ©, 
A-« J B(M 0> AT([0,At]) 9 )^ J B(M ( J([0,T])«) 

By choosing T = T'A -1 and using the identity 

M (S-iT<([°> A ~ lr ]) ~ aW *~ lM o%>([°, T '}), we in- 
clude that 

E(M o ,T'([0,Xt\y) ~ A f « 

A — >0 

This achieves the proof. 

We can therefore see, that, as far as asymptotic mul- 
tifractality is concerned, the pertinent parameter is the 
small scale behavior of the function f(l) or equivalently 
the small time behavior of pi—o(t). As pointed out pre- 
viously, in the case of a Gaussian field u>(t) (i.e. the 
infinitely divisible law has only a Gaussian component), 
p(t) is nothing but the covariance of the process. The 
previous discussion leads thus to the conclusion that 
non trivial limit multifractal measures arise only in the 
marginal situation when the correlation function of the 
logarithm of the fluctuations decreases as a logarithmic 
function. 

G. An alternative discrete time construction for 
MRM 

In the case J x~ 2 v(dx) < oo (e.g. the Levy mea- 
sure has no mass in an interval around x = 0) a real- 
ization of the measure P(dt, dl) is made of dirac func- 
tions distributed in the iS + half-plane. Thus the process 
M;(([0, t]) = f e u '^dt is a jump process that can be 
simulated with no approximation. 
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However, if x~ 2 v(dx) has a non finite integral (e.g., it 
has a Gaussian component), this is no longer the case. 
Thus, one has to build another sequence of stochastic 
measures Mi(dt) that converges in law towards M(dt) 
and that can be seen as a discretized version of M(dt). 
We will see, in following sections, that such a discrete 
time approach is also interesting for multifractal stochas- 
tic processes construction. 

We choose Mi(dt) to be uniform on each interval 
[kl,(k + l)l[, Vfc £ N and with density e ""W. Thus, 
for any t > such that t — pi with p £W , one gets 



p-i 



M J ([0,t]) = E e 



(63) 



k=0 



In the same way as for the measure Mt(dt), one can prove 
f|, within the framework of positive martingales, that, 
almost surely, M l=2 -™{dt) converges towards a well de- 
fined limit measure when n — > +oo (i.e., Z — » + ). More- 
over, in the same way as in section III B, if we suppose 
■0(2) < 1, then one can show (the proof is very similar as 
the one in Appendix A) that, in the mean square sense, 



lim M l=2 - n {dt) = M(dt). 



(64) 



As long as i/j(2) < 1, this construction gives therefore a 
way of generating a measure which is arbitrary close (by 
choosing I small enough) to the limit measure M(dt). 



IV. LOG-INFINITELY DIVISIBLE 
MULTIFRACTAL RANDOM WALKS 

In this section we build and study a class of multi- 
fractal stochastic processes that are no longer, as before, 
strictly increasing processes (measures). They can be ba- 
sically built in two different ways: (i) By subordinating 
a fractional Brownian motion (fBm) with the previously 
defined MRM M (t) or (ii) by a stochastic integration of 
a MRM against a fractional Gaussian noise (fGn). As 
we will see, most of the statistical properties of the so- 
obtained random processes are directly inherited from 
those of the associated MRM. 



where dW(t) is a Gaussian white noise independent of 
u)[ . The MRW is then defined as the limit of X\ (t) when 
1^0+ : 



X(t) = lim XAt). (66) 

1^0+ 



One can easily prove that for fixed t and I, one has 



Xit) 1 ^ m {t)e (67) 



where e is a standardized normal random variable inde- 
pendent of 07 it) which is itself nothing but the associated 
MRM as defined previously: 



s Mv)du = Mi(t) 



(68) 



Let us note that the (non decreasing) increments of of (t) 
is referred to (in the field of mathematical finance |2^, 
p0| ) as the stochastic volatility. 

Using the same kind of arguments on finite dimensional 
laws, one can also prove that the finite dimensional laws 
of the process Xi (t) converge to those of the subordinated 
process B(M(t)). Actually, one can show || that, as 
long as ip'(i) < 1 ( i.e., Be > 0, 0(1 + e) < 1 which is, 
as mentioned in section II B, the condition for the limit 



measure M{dt) to be non degenerated), on has 

X(t) = lim X l (t) l = lim B(ai(t)) l = B(M(t)). (69) 

1^0+ 1^0+ 

The so-obtained MRW can be thus understood as a 
Brownian motion in a "multifractal time" M(t). The 
subordination of a Brownian process with a non decreas- 
ing process has been introduced by Mandelbrot and Tay- 
lor |Jl| and is the subject of an extensive literature in 
mathematical finance. Multifractal subordinators have 
been considered by Mandelbrot and co-workers |52| and 
widely used to build multifractal processes from multi- 
fractal measures (see below). In a forthcoming section 
we will see that multifractal subordination and stochastic 
integration do not lead to the same processes when one 
considers long range correlated Gaussian noises (fGn). 



A. MRW with uncorrelated increments 

In the same spirit as the log-normal MRW construc- 
tions in Refs. |l|, 0] (see also [^?J), we use a stochas- 
tic integration of e^'- 1 against the (independent) Wiener 
measure dW(t). 

1. Definition 

bmce one can consider the process 

Xi(t) = f e^ l{u) dW{u) (65) 
Jo 



2. Expression of the moments and multifractal properties 

Thanks to Eq. ( |69|) (assuming < 1), one gets the 

expression of the absolute moments of X(t) (or X(to + 
t)-X(t )): 

Vg, E(\X(t)\i) = E{\e\*)E(M{[V,t])^) 

= °* r( i) a W M([0 ' t]) ) ' 

where the first factor comes the order q moment of a 
centered Gaussian variable of variance a 2 . 
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If M(dt) is an exact multifractal stationary random 
measure, then, X(t) obeys the exact multifractal scaling 
equation: 

E (\X(t)\«) = a" t{ \ i > K q/2 t^ (70) 



where K q is defined in Eq. ( f43| ) and 

C q = q/2 - TP(q/2) . 



(71) 



Using Eqs. @, and @, one deduces that, for 
q> 2, 

S(|X(t)|«)<+oo^C«>l, (72) 

and conversely, 

C <? <l^ J B(|X(t)|5)<oo. (73) 

Moreover, let us note that if M (dt) verifies only an 
asymptotic scaling, so does the MRW process X(t). 



3. An alternative discrete time construction for MRW 
processes 



As in Refs. |jj |2| or in section IIIG, one can also try 
to build an MRW process using a discrete approach. A 
discrete construction can be useful for numerical simu- 
lations. Let us, for instance, choose l n = 2~ n . And 
let ei[k] = dW{u) be a discrete Gaussian white 

noise. We define the piece- wise constant process Xi n (t) 
as (t = pl n ): 



*!„(*) 



p-l 

y^ e |o) !n (fei„) 

k=0 



(74) 



The MRW Xi n (t) can be rewritten as 



p-i 



fc=0 



where Mi n is defined by (63). One then deduces easily the 
convergence of Xi n from the convergence of Mi n . Thus, 
for instance, one can prove || that, as long as ip{2) < 1, 
one has 



lim Xr (t) = B(M(t)). 



(76) 



B. MRW processes with long-range correlations 

1. Definitions 

In order to construct long-range correlated MRW, it 
is natural to replace the Wiener noise (resp. Brownian 



motion) in previous construction by a fractional Gaus- 
sian noise (resp. fractional Brownian motion). A fBm, 
Bn(t) is a continuous, self-similar, zero- mean Gaussian 
process which covariance reads (see e.g. for a precise 
definition and properties): 



E(B H (t)B H ( S )) = Y ( S 



2H + t 2H -\t-s\ 2H ) 



(77) 



where < H < 1 is often called the Hurst parameter. 
Standard Brownian motion corresponds to H = 1/2. 

The simplest approach to construct a long-range cor- 
related MRW follows the idea of Mandelbrot that simply 
consists in subordinating a fractional Brownian motion 
of index H with the MRM M(t), i.e., 



X° H (t) = B H [M(t)\ 



(78) 



An alternative would consist in buidling a stochastic 
integral against a fGn dWfl-(i) 



(79) 



and considering some appropriate limit I — * 0. How- 
ever stochastic integrals against fGn cannot be defined as 
easily as for the white Gaussian noise and the proposed 
constructions require the complex machinery of Malliavin 
calculus or Wick products [p3|, EjJ) . One simple way to de- 
fine the previous integral could be to see it as the limit 
of a Ricmann sum: 



t/At 

e^ (t) dW H {t) = lim V e"' 

„ At^O ^— ' 



(kAt) e H Mk] (80) 



where e H ,At(k) = B H (kAt) - B H ((k - l)At). We have 
not proved yet that this is a mathematically sound defini- 
tion. However, if one assumes that (|79| ) makes sense, one 
can address the question of the existence of the limit pro- 
cess lim;^o H (t). In Appendix B, we provide heuristic 
arguments for mean square convergence. We obtain a 
condition: 



H > l/2 + V(2)/2 



(81) 



where ip(p) is the cumulant generating function associ- 
ated with oj. 



2. Multifractal properties 

In the case of the subordinated version Xfj(t) of the 
MRW, the scaling properties can be directly deduced 
by the self-similarity of B H (t) Q. Since B H (t + r) - 
B H (t) = law t h (B H {t+l)-B H {t j), and M(t) is indepen- 
dent of B H , one has X H (M(t)] \= law M(t) H B H (l)- The 
scaling of the absolute moments of the increments of Xh 
is therefore: 

E[\X H (t)\«] = E[M(ty H ]E[\B(l)n 
= K qH G q t^ 
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with 



C =qH- 4>{qH) 



(82) 



For the second version, X l H , the scaling of the moments 
is determined using the scale-invariance of the process 
u>i(t) and the self-similarity of the fGn dW#(i). Using 
the same method as for the measure M (t), one obtains: 



E[\X* H {t)\i] 



Mgt C » 



with 



C = qH- i>(q) 



(83) 



(84) 



We can see that the multifractal spectra of X S H and X l H 
are different: they do not correspond to identical pro- 
cesses as it was the case for the uncorrelated construc- 
tion. Notice that the existence criterion ( |Sl| ) , can be sim- 
ply rewritten as Q > 1- A ccording to the considerations 
developed in section III D| , this condition ensures the ex- 



istence of the second order moment of X l H . Wether the 
class of processes X l H (t) can be extended, in some weak 
probabilistic sense, to values of H < 1/2 + -0(2)/2 is still 
an open problem (such processes would have an infinite 
variance). The condition of finite variance for Xfj(t) is 
less restrictive since it comes to the condition K 2 h < oo, 
where K 2 h is defined in Eq. (|43|). For H < 1/2 such a 
moment is always finite. 



C. A remark on subordination 

Let us remark that, in some sense, the subordination 
by a MRM M(dt) can be iterated. Indeed, if Mi([0,t]) 
and M 2 ([0,i]) are two independent MRM, the subordi- 
nated measure 



M([t 1 ,t 2 ]) = M 1 (M 2 ([0,t 1 ]),M 2 ([0,t 2 ])]) 



(85) 



is well defined. Using the cascade equation (|40|), we de- 
duce that 

M([0,Xi\) = law Mi (o,Wf } M 2 ([0,t]) (86) 
= law w£] a) M([0,t]) (87) 

( 1 21 

where ' are the (independent) log-infinitely divisi- 
ble weights associated with the MRM Mi i2 . The second 

(2) 

equality is valid only if W x < 1, i.e., when the Levy 
measure associated with M 2 is concentrated on ]0, +oo]. 
By computing the moment of order q of both sides of 
the equality, we see that the mulifractal spectrum of the 
subordinated measure reads: 



c, 



A2) 

S>9 



^(q-^Xq)) (88) 



where (resp. ip^(q)), i — 1,2, is the spectrum asso- 
ciated with Mi. The equation (87) corresponds to a "ran- 
domization" of the rescaling factor A that parametrizes 



the log-infinitely divisible law of ln(W_j[ ). It is easy to 
prove || that the law of \n(W {1) (2} ) remains infinetely di- 
visible. The class of log-infinitely divisible MRM is there- 
fore closed under subordination. The family of subordi- 
nated spectra ( pq ) is thus included in the family of log- 
infinitely divisible spectra and the operation ( p5[ ) does 
not allow us to build new MRM with exact scale invari- 
ance properties. 



V. NUMERICAL SIMULATIONS 
A. Principles 

In order to generate realizations of Mi n (dt) (defined 
by (|63[)), one needs to be able to generate realizations 
of u>i n (kl n ). However, following the definition of uoi 
one would need realizations of the 2d random measure 
P(dt,dl) for I > l n . In the case of a compound pois- 
son process, the process M;([0, t]) = J* e^'^dt is a jump 
process that can be synthetized easily. 

In the general case, we need to find a set of disjoint "el- 
ementary" domains of the half-plane S + such that, for 
any k, there exists a subset of this set such that A\ n (kl n ) 
can be expressed as a the union over the elementary do- 
mains of this subset. Since, at fixed l n , the boundaries 
of the domains Ai n (kl n ) (k e 7L) define a tiling of S + , it 
is natural to consider the elementary cells of this tiling. 
Each cell is the intersection between left and right strips 
limited by left and right boundaries of conical domains: 
Let us define the cell Bi(t,t') (with t < t') as 

8^,0 = (Ai(t)\Ai(t - I)) n (Ai(t')\Mt' + 0) (89) 

Then, by definition the cells {Bi n {kl ni k'l n )}k<k' are dis- 
joint domains and form a partition of the subspace (of 
S+) {(t,l) ES+, l> l n }. Moreover 

A n (M n )= (J (J B ln (il n ,jl n ). (90) 

— oo <i<k /c<j < + oo 

On the other hand, for a fixed s>l ni one has 

(u,s) e Ai(t)\Ai(t-l) «. t-l + f(s)/2 <u< t+f(s)/2, 

and 

(u,s) e Ai(t + l)\Ai(t) & t-f(s)/2 <u< t+l-f(s)/2. 

Thus setting Yij = P(Bi n (il n , jl n )), straightforward 
computations lead to the following representation of dis- 
crete process LUi n (kl n ): 



k +OG 

i— — oo j—k 



(91) 



where {^j jij are independent infinitely divisible ran- 
dom variable which satisfies 



E (e^- 



(92) 
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FIG. 3: MRW signals sampled at rate 1. (a) log-Normal 
MRW with T = 512 and A 2 = 0.05. (b) log-Poisson sample 
with T = 512, 7 = 4 and 7ln(<5) 2 = 0.05. In both signals, the 
cut-off I as been fixed to 1/8 for the numerical synthesis. 



with 



f ds/s 2 [ dtH iij<s {t), 

Js>L, J 



(93) 



where Hij tS (t) is the indicator function 

Hi,jAt)=\a i , j (s),b i , j (s)](t), (94) 

with 

a itj (s) = max {(i - l)l n + f(s)/2Jl n - f(s)/2) , (95) 
and 

b itj (s) = min {il n + f{s)/2, (j + l)l n - /(a)/2) . (96) 

Let us note that if the function / is bounded (which is 
the case if we are under the hypothesis (p2|), i.e., in the 
case of "exact" scaling) , the number of terms in Eq. ( pl| ) 
is finite. 

In the Gaussian case, the situation is simpler because 
mi is completely characterized by its covariance function, 
Pi(t). In that case, the 2D synthesis problem can be 
easily transposed as a ID filtering procedure by finding 
a filter (f>i(t) such that 

b*4n=Pi, (97) 

where * stands for the convolution product. The process 



U - 



FIG. 4: £q spectrum estimation for log-Normal and log- 
Poisson MRW's. The exponents have been directly estimated 
from linear regression of increment order q absolute moments 
in doubly logarithmic representations. These moments have 
been estimated using a statistical sample of 256 signals of 64 
integral scales. The symbols (•) correspond to log-Poisson 
estimates while (o) correspond to log-Normal estimates. The 
solid line represents the log-Poisson analytical spectrum ( 55 ) 
and the dashed line to log-Normal analytical spectrum (51) 
rescaled using Eq. ([n]). The parameters of both processes 
are those of Fig. 2 and have been chosen so that the so called 
intermittency parameter —(""(0), is the same (0.05). 



B. Numerical examples 

In Fig. 3 are shown two samples of MRW which are re- 
spectively log-Normal and log-Poisson (see Eqs. ( ^l| , p2[ )). 
In both cases we have chosen T = 512 sample units, 
if = 1/2 (e is a Gaussian white noise). For the log- 
Normal process, A 2 = 0.05 while 7 = 4 and S = e -V°- 05 /7 
for the log-Poisson process. The log-Poisson process has 
been synthetized using Eq. (|9l]) while a simple filtering 
method was used for the log-Normal process. Eq. ( |97| ) 
was solved numerically in the Fourier domain. In Fig. 4 
are plotted the ( q functions estimated for both processes. 
These functions have been obtained from the scaling of 
the moments estimated using 256 MRW trials of 64 inte- 
gral scales long. The ( q values for negative q have been 
obtained using the so-called WTMM method that is a 
wavelet based method introduced to study multifractal 
functions |3f| [36| |37j. The superimposed analytical for- 
mulae obtained with Eqs. (fn]), (|l|) and (||), fit very 
well with statitical estimates except a large negative q 
values for the log-Normal case. This can be explained as 
a finite statistics effect. 



w,'(t) = / 4>i(t-t')w{dt') 



(where W is here a ID Wiener noise) will be thus identical 
to the normal process u)iCt) as defined in Eq. (|2TJ) . In the 
original study of Refs. |ll ^ the MRW has been defined 
along this method. 



VI. CONNECTED APPROACHES 

In this section we review some specific results concern- 
ing respectively log-normal and log-Poisson compound 
MRW. 
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A. Log-normal MRW 

The log-normal MRW has been originally defined in 
Refs. f|, §]■ It corresponds to the simplest situation 
when the Levy measure has only a Gaussian component. 
In that case, the Gaussian process u>i(t) can be directly 
constructed from a ID white noise, without any refer- 
ence to 2D conical domains. This model is interesting 
because its multifractal properties are described by only 
two parameters, the integral scale T and the so-called 
intermittency parameter A 2 . Some simple estimators of 
these quantities have been proposed in jjj. Moreover, 
many exact analytical expressions can be obtained and 
notably the value of the prefactor K q in Eq. (^) . In |l0| , 
it is shown that this prefactor can be written as a Selberg 
integral Elt] . Its analytical expression reads: 



9-1 



K„ 



k=0 



r(i-2A 2 fc) 2 r(i-2A 2 (fc + i)) 

r(2 - 2X 2 (q + k - l))r(l - 2A 2 ) 



(99) 



It is easy to check that K q is defined only if q < q* = 
2/ A 2 . We recover the finite moment condition (49), 
q < q* with Q qt — 1. Notice that, as emphasized by Frisch 
p2[ , the existence of infinite moments can be a drawback 
of a log-normal multifractal as a model for experimental 
situations like turbulence. However, for a typical value 
A 2 = 5.10 -2 , = 40, so a log-normal approximation 
can be very good in a range of q values far beyond the 
limit associated with the finite size of experimental sam- 
ples. Let us finally mention that the log-normal MRW 
can be naturally generalized to a "multivariate multi- 
fractal model" which is a multifractal vector of processes 
characterized by an intermittency matrix Ay Q. This 
notion of "joint mult ifr act ality" can be very interesting 
in many applications. 



B. Multifractal products of cylindrical pulses 

In reference (|], Mandelbrot and Barral introduced a 
positive multifractal random measure using products of 
positive random variables associated with Poisson points 
within 2D conical domains. Their construction is a par- 
ticular case of MRM. It actually reduces to the case where 
v{dx) is a Levy measure satisfying J v(dx)x~ 2 < oo 
(see section III E| ) and where the set Ai(t) is delimited 
by the function /(/) as defined by Eq. (jl7). Let us 
note that, since they have not considered the full do- 
main Ai(t) associated with the fun ction f(l) as defined 
by Eq. (|22j), as explained in section III F , this construc- 
tion performs asymptotic scaling (ji|) but not exact scal- 
ing ([!]). However, these authors did not study the scaling 
properties of the random measures. They rather focused 
on the pathwise regularity properties. More precisely, 
they proved the validity of the so-called "multifractal 
formalism" (see e.g. @ |§ |7j [3§ ||) that relates 
the function £ 9 to the singularity spectrum D(h) associ- 
ated with (almost) all realisations of the process. For a 



given path of the increasing process associated with an 
MRM, D{h) is defined as the Hausdorff dimension of the 
set of "iso-regularity" points, i.e., the points where the 
(Holder) regularity is h. Barral and Mandelbrot proved 
that D(h) and Q q are related by a Legendre transforma- 
tion. Since we proved that £ g is the scaling exponent 
of MRM moments, it follows that one can estimate the 
singularity spectrum of the MRM paths in the case of 
log-Poisson compound statistics. It should be interesting 
to extend the Mandelbrot-Barral theorem to the general 
log-infinitely divisible MRM and MRW paths. 



VII. CONCLUSION AND PROSPECTS 
A. Summary and open questions 

In this paper we have constructed a class of station- 
ary continuous time stochastic measures and random pro- 
cesses that have exact or asymptotic multifractal scaling 
properties in the sense of Eqs. (Q) and (H). We have 
shown how stochastic integration of an infinitely divisi- 
ble noise over cone-like domains, as originally proposed 
in Q, naturally arises when one wants to "interpolate" 
discrete multiplicative cascades over a continuous range 
of scales within a construction that is invariant by time 
translations. The exponential of these stochastic inte- 
grals (e" !< - t - ) ) can be interpreted as a "continuous prod- 
uct" from coarse to fine scales and thus as the continu- 
ous extension of the multiplicative rule involved in the 
definition of discete cascades. We have shown that the 
probalility density functions associated with MRM and 
MRW processes can have, like discrete cascades, fat tails 
with arbitrary large exponents. However, unlike their dis- 
crete analog, our "continuous cascades" have stationary 
fluctuations, do not involve any particular scale ratio and 
can be defined in a causal way. This "sequential" formu- 
lation, as opposed to the classical "top to bottom" defi- 
nition of multifractals, can be very interesting for mod- 
elling dynamical processes (see the next section) . Let us 
note that we focused in this study on ID processes but 
our construction can easily be extended to higher dimen- 
sions. 

It is well known (see Refs. |l], |j) that the multiscaling 
(|l|) or (^J) with a non-linear convex Q q function cannot ex- 
tend over an unbounded range of scales and there neces- 
sarily exists an "integral scale" T above which the scaling 
of the moments changes. The existence of such an inte- 
gral scale can be found in the g eneral shape of function 
f(l) as discussed in section IIIF: One must have /(/) ~ I 
when I — > and f(l) = o(l) when I — > +oo, so the scale 
T is a scale that separates these two asymptotic regimes. 
It is remarkable that there exists a particular expression 
for the function f(l) (Eq. (|22])) for which the moments 
in the multifractal regime I < T satisfy an exact scaling. 
The existence of processes with such properties was not a 
priori obvious. In the same section, we have shown that 
f(l) ~ I is a necessary and sufficient condition for the 
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existence of a limit multifractal object. From a funda- 
mental point of view, one important question concerns 
the unicity of our construction: Is any process satisfy- 
ing (Q) can be represented within the framework we have 
introduced ? 

It remains many open mathematical problems related 
to the processes which we introduced in this paper. Some 
of them have been already mentionned, notably the ques- 
tions re lated to the construction of sto chast ic integrals in 
section [VB. As discussed in section VI B, it should be 
interesting to generalize the results of Ref. || in order 
to link scaling properties and pathwise regularity within 
a multifractal formalism. Another interesting problem 
concerns the study of limit probability distributions asso- 
ciated with MRM for which fery few features are known. 
Like infinitely divisible laws, they appear to be related 
to some semi-group structure. Finally, one can wonder if 
log-infinitely processes we have defined are not the nat- 
ural candidates to be described within the framework of 
"Markovian continuous cascades" as introduced in Refs. 

MM. 



B. Possible applications 



analogy between turbulence and finance |5l], |5^, |5J, g5| . 
Logarithmic decaying correlations and "1/f" power spec- 
trum have been directly observed for various time series 
f34| |56| , so it is reasonable to think that MRW models are 
well suited for modeling financial time series W . The ver- 
satility of infinitely divisible MRW is very interesting to 
account for various stylized facts of financial times series 
such that the multiscaling, the power-law tail behavior of 
return pdf and therefore such models can be very helpful 
for financial engineering and risk management. 

Among all the remaining disciplines where multifrac- 
tality has been observed, one can mention the study of 
network traffic j5?], |5^, p9|j , geophysics and climatology 
p0| . |6l| , |62| , biomedical engineering |33), the modelling 
of natural images [|64[ , . . . On a more theoretical ground, 
"continuous branching trees" and log-correlated random 
processes have also been considered in the physics of dis- 
ordered systems Q . The study of free energy density 
in the thermynamic limit in presence of log-correlated 
disorder (or equivalently on topologies involving random 
disordered trees) raise questions very similar to the study 
of limit MRM addressed in this paper. One can hope that 
pushing forward this analogy in a promising prospect to 
get significant results in both fields. 



One of the main issues of the present work was to con- 
struct a wide family of multifractal processes (or mea- 
sures) that are likely to be pertinent models in many 
fields were multiscaling laws are observed. Naturally, the 
first application of which one can think, is fully devel- 
oped turbulence. Turbulence and multifractals share a 
long history and we refer the reader to Ref. |^2j for a 
review on the "intermittent" nature of turbulent fields. 
Recently, new aspects of turbulence were studied by con- 
sidering fluid dynamics from a Lagrangian point of view. 
This was possible because two groups developed new ex- 
perimental devices based on a fast imaging system Q 
or ultrasound techniques |^6|, [l7j allowing for a direct 
measurement of the velocity of a single tracer in a turbu- 
lent flow. In a recent works, Pinton and his collaborators 
Pq, [ItJ studied the intermittency of Lagrangian trajecto- 
ries and related it to the slow (logarithmic) decay of the 
particle acceleration correlations j|7|, , very much like 
for a MRW model. In other words, these authors found 
that the turbulent Lagrangian dynamics is very well de- 
scribed by an equation of Langevin type with a driving 
force amplitude similar to e W( W involved in the MRW def- 
inition. The understanding of the physical origin of such 
dynamical correlations and the link between Lagrangian 
and Eulcrian statistics is a very promising path towards 
the explanation of the intermittency phenomenon in fully 
developed turbulence. 

Besides turbulence, "econophysics" |5(| |||, [32| is an 
emerging field where fractal and multifractal concepts 
have proven to be fruitful. Indeed, many recent studies 
brought empirical evidences for the multifractal nature 
of the fluctuations of financial markets (see jl| and ref- 
erence therein). Some physicists raised an interesting 
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APPENDIX A: MEAN SQUARE 
CONVERGENCE OF A MRM 



Let Mi(dt) defined as in Eq. (p3|). In this section we 
prove that, assuming ip(2) < 1, one has 



M,([0,t])^M([0,t]) 



Let us define 



i? M ,( T ) = E Lf->iV)+»A«+r)\ 



(Al) 



(A2) 



In order to prove (Al), let us first show that, if I' < I: 

Rl,v(r)=Ri, l (T) = e"^^ (A3) 



with pi(r) as defined as in ( p0[ ) and (|3l|). 

The first equality in Eqs. (A3) comes directly from 
the assumption — while the second equality is a 
particular case of the identity (|33]), where q = 2, p\ = 
P2 — —i and t\ — t% = r. 

Let us show, that, V e, 3 Iq, V 1,1' < Iq, 
E [{Mi(t) - M r (i)) 2 ] < e. Let us suppose that I' < I. 
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Then, 



of the process 



C(l,l',t) = E [(M,([0,t]) - M r ([0,t])) 2 ] 
= B[M, 2 ([0,t])] +f?[M,?([0,i])] - 2S[M«([0,t])M I /([0,i])] 
(e ( e Mu)+"i(.v)\ +E ( e "i>M+"i>(v))) dudv 



X t (t) = / dW„(t)e"^ 
Jo 



(Bl) 



'o J o 
Jo Jo 



Thus, thanks to Eq. (|A3|), we get, after a little algebra 



C{l,l',t) < Dt (R v>v {u)-Ri,i(u))du, 



< Dt f R v .v{u)du, 
Jo 



< Dt [ R v .v{u)du + Dt ( R v . v {u)du 
Jo Jv 

,/l-V>(2) 



< Etr 



+ Et (l 1 -*® - Z' 1 -^) 



where D and £ are positive constants. Since ^(2) < 1, 
we see that M;([0, t]) is a Cauchy sequence and thus con- 
verges in mean square sense. 

The previous computations also prove the conver- 
gence of all finite dimensional vector {Mi[ti,t\ + 
Ti], . . . , Mi[t n , t n + T n ]}. In order to prove the existence 
of the limit M([0, t]) as a stochastic process, one further 
needs a "tightness" condition. Such a condition can be 
obtained along the same line as for previous computation. 
Indeed, provided ip{2) < 1, one can bound the order 2 
moment of Mj([0, t]), Vt, I: 

E(M l ([0,t] 2 ) < Ct 2 -^ , 

where C does not depend on I. This achieves the proof. 



APPENDIX B: MEAN SQUARE CONVERGENCE 
OF A MRW WITH A FRACTIONAL GAUSSIAN 
NOISE e H 

In order to simplify the proof and to avoid technical 
complications, let us show the mean square convergence 



where dWu is a continuous fGn which covariance is (H ^ 
1/2) 



lH ( T ) = a 2 H(2H-l)r 



2H-2 



(B2) 



Rigorously speaking the previous integral is not well 
defined but the proof of the convergence of the discrete 
version (p0|) is very similar. Let us now 



(B3) 



Xl (t)^X(t) 

provided H satisfies: 

H > l/2 + V(2)/2 



(B4) 



We proceed alo ng the same line as in Appendix A. In 
order to prove (B3), we have to prove, that, V 1,1' < I, 
E[{Xi(t)-Xy{^f] -j_*o0. Let 

E = E[(X l (t)-X l ,(t)) 2 ] . (B5) 
Thanks to Eq. (|A3|), we have, 
E = E[X?(t)]+E[X?,(t)]-2E[X l (t)X l ,(t)} 

= Jo Jo ( E ( eMu)+Mv) ) + E (eP'< u '> +u *'Wy) j H (\u-v\)dud 
2 ! [ E { e "iW+"i>( v A lH (\ u - v\)dudv 



JO 
t r t 



Jo 



(Ri'i'Qu - v\) - Ri,i{\u - v\)) jh(\u - v\)dudv 



The last integral behavior in the limit I — > can be easily 
evaluated. After some simple algebra we get: 



E[{x x {t)-x v {t)Y] =o(i 2H -^- 1 



Thus if condition (B4) is satisfied, i.e. 2H-tp(2)-l > 0, 



Xiit) is a Cauchy sequence and thus converges in mean 
square sense. 
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